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Abstract: A construction is presented that allows to produce subspace codes of long length 
using subspace codes of shorter length in combination with a rank metric code. The subspace 
distance of the resulting code, called linkage code, is as good as the minimum subspace distance 
of the constituent codes. As a special application, the construction of the best known partial 
spreads is reproduced. Finally, for a special case of linkage, a decoding algorithm is presented 
which amounts to decoding with respect to the smaller constituent codes and which can be 
parallelized. 
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1 Introduction 

In [22] Koetter and Kschischang developed an approach to random network coding where the 
encoded information is represented as subspaces of a given ambient space. This accounts for the 
unknown network topology by assuming that any linear combination of packets may occur at 
the nodes of the network. 

This approach has led to the area of subspace codes and specifically to intensive research 
efforts on constructions of subspace codes with large subspace distance [23l|2ll|25l[ini[ia[9lll6l 
ttSKHKIIllMlEHKlT]. Most of the research focuses on constant-dimension codes (CDC’s), that 
is, codes where all subspaces have the same dimension. 

One direction for constructing CDC’s is based on so-called cyclic orbit codes [23l |9l [241 
EIKITI, which are orbits of a subspace in the Fg-vector space F^n under the natural action of 
F*. While the resulting codes have very beneficial algebraic structure, they do not have large 
cardinality in general. Taking unions of such codes leads to cyclic subspace codes which still 
have nice structure, but it remains an open problem how to take unions of cyclic orbit codes 
without decreasing the distance. 

A second major research direction is based on rank-metric codes as introduced and studied 
earlier by Delsarte [6] and Gabidulin m- Lifting rank-metric codes [26] is a very simple con¬ 
struction which results in subspace codes where the reduced row echelon form of each subspace 
has its identity matrix in the leftmost position. While these codes are asymptotically good [21], 
they can still be improved upon. Through a careful study of general reduced row echelon forms, 
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Etzion and Silberstein [10] could enlarge lifted rank-metric codes. This is known as the echelon- 
Ferrers construction (or multilevel construction) and has led to the best codes known at that 
time. In the same paper it is also demonstrated how to decode these codes. In m the same 
authors could enlarge their codes even further by making use of the pending dot construction 
initiated by Trautmann/Rosenthal [30]. A different way of enlarging lifted rank-metric codes 
has been presented by Skachek in [27|, where for a specific case of the resulting codes also a 
decoding algorithm is developed. 

A slightly different approach is taken in [12( Sec. Ill] by Etzion/Vardy and [l9| by Gorla/Ra- 
vagnani. Therein, the authors focus on CDC’s with the best theoretical distance. This distance 
is achieved when any two distinct subspaces of the code intersect trivially. In |12[ [19] the authors 
present constructions of such codes with large cardinality (in fact, the same cardinality). In finite 
projective geometry, collections of pairwise trivially intersecting subspaces are known as partial 
spreads, or as spreads if the subspaces also cover the entire ambient space. The cardinalities 
obtained in naiHi have also been achieved via constructions based on finite geometry [D El- 

In this paper we will present a surprisingly simple construction of subspace codes of long 
length by linking matrix representations of subspace codes of shorter length and making use of 
a rank-metric code. The distance of the resulting subspace code, called a linkage code, is the 
minimum of the distances of the involved subspace codes and the related lifted rank-metric code. 

By way of an example we will compare our construction to the echelon-Ferrers construc¬ 
tion m and its improvement m as well as to the best codes found by sophisticated computer 
search. We will see that our codes improve upon those constructed in mm and for lengths 
where no computer search has been conducted, it leads to the best codes known so far. 

In Section [3| we show that specific instances of linkage reproduce the constructions of partial 
spreads given in imiini, and thus linkage may be considered as a generalization of the latter. 
In fact, those constructions, which on first sight appear to be quite different, can be nicely 
presented in a unified way with the aid of linkage. Using an optimal binary partial spread of 
dimension 3 and length 8, which was found only recently by El-Zanati et al. |8|, we can improve 
upon the constructions in [121 [19] for the binary case in dimension 3 and when the length has 
remainder 2 modulo 3. This reproduces results derived in [8]. 

Finally, in Section 0] we show that for a specific choice of the seed codes, decoding of linkage 
codes can be reduced to decoding of the seed codes. Making use of Gabidulin codes and lifted 
Gabidulin codes, for which an efficient decoding algorithm has been derived by Silva et al. [26j . 
we obtain linkage codes that can be decoded by an efficient parallelizable algorithm. 

2 A Linkage Construction for Subspace Codes 

Let us first recall some basic facts about subspace codes and rank-metric codes. Throughout 
we fix a finite field F = Fg. A subspace code of length n is simply a non-empty collection of 
subspaces in F*^. The code is called a constant-dimension code if all subspaces have the same 
dimension. The subspaee distanee of a subspace code C is defined as ds(C) := min{ds(V, W) | 
V, W € C, V / W}, where the distance between two subspaces is 

ds(V,>V) :=dimV + dim>V-2dim(Vn>V). 

It is a metric on the space of all subspaces, see [221 Lem. 1]. If dimV = dimW = k, then 
ds(V, W) = 2(dim(V -|- W) — k). As a consequence, if C is a constant-dimension code of 
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dimension k then 


ds(C) < min{2/c, 2(n — k)}. (2-1) 

If the code C consist of a single subspace of dimension k, we define ds(C) := min{2A:, 2{n — k)}. A 
constant-dimension code of length n, dimension k, cardinality N will be called an (n, N, k)-code, 
and it is a (n, N, k, d)-code if its subspace distance is d. 

A k X m rank-metric code is a non-empty subset of endowed with the rank metric 

dR(A, i?) := rk(A — B) (which is indeed a metric, see |13]1. The rank distance of a rank-metric 
code C is defined in the usual way as dR(C) := min{rk(A — B) \ A, B ^ C, A ^ B}. If C consists 
of single matrix, we define dR(C) := min{/c, m}. It is well known (see [U Thms. 5.4, 6.3] and |13l 
p. 2]) that if m > k and C is a rank-metric code in with rank distance d, then 

\C\ < qm{k-d+i)_ ^2.2) 

Moreover, there exist rank-metric codes of distance d and size qi^ik-d+i) ^ such codes can 
even be constructed as linear subspaces of F^^”*. They are called MRD codes. The best known 
class of linear MRD codes are the Gabidulin codes, derived by Gabidulin in [13] . Just recently, 
other constructions of MRD codes were found by de la Cruz et al. |5] and Hernandez/Sison |20j . 
However, as opposed to Gabidulin codes so far no decoding algorithm is known for the latter 
codes. 

If C C F^^™ is an MRD code, then the subspace code 

C = {im(/fc I A) I A € C} 

is called a lifted MRD code. Here the notation im(M) stands for the row space of the matrix M 
and Ik denotes the k x /c-identity matrix. If dR(C) = d, then ds(C) = 2d, see |26l Prop. 4], and 
therefore C is a (A; -|- m, qM^-d+i) ^ 2d)-code. 

The following specific class of MRD codes will be used in the next section when studying 
partial spreads. 

Remark 2.1. There is a simple construction of MRD codes in F^^"* of rank distance k. Let 
W G be any matrix of rank k and M G GLm(F) be the companion matrix of a primitive 

polynomial in F[a;] of degree m. Define 

C := {WM^ I Z = 0,..., - 2} U {0} C 

Since F[M] = F^m and \M\ = — 1, the code C is a linear rank-metric code of size and 

rank distance k. Hence C is an MRD code. In fact, it can be shown that C is a Gabidulin code. 

In order to present our linkage construction we need to work with matrix representations of 
subspaces. The following terminology will be convenient. 

Definition 2.2. A set of matrices A4 C F^^” is called SC-representing if rk(M) = k for 
all M G A4 and im(M) ^ im(M') for all M ^ M'. The induced constant-dimension code 
{im(M) I M G M.} is denoted by C(M). 

For example, the rank-metric code C in Remark 12.11 forms an SG-representing set of a cyclic 
orbit code in the sense of ISlEHIttTI. Any set of full row rank matrices in reduced row echelon 
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form is an SC-representing set. In general, an SC-representing set is simply a subset of orbit 
representatives of the action of GLfc(F) on via left multiplication. 

The linkage construction in the following theorem links subspace codes with the aid of a 
rank-metric code and results in a subspace code of longer length without compromising the 
distance. It makes use of representing matrices. The theorem generalizes a former construction 
in [171 Thm. 5.1]. 

Theorem 2.3. For i = 1,2 let M.i C he SC-representing sets of cardinality Ni. Thus Ci = 

C{M.i) is an {ui, Ni, k)-code. Let ds(Ci) = di. Furthermore, let Cr C F^^”^ ^ linear rank- 

metric code with rank distance dR(CH) = dn and cardinality \Cji\ = Nr. Define the subspace 
code C of length n := ni + n 2 as C := Ci U C 2 U C 3 , where 

Cl = {im(t/ I Ofcxna) I U € Adi}, 

C 2 = {im(0fcxni \U)\U e M 2 }, 

C 3 = {im{U\M)\U gMi, M G Ci?\{0}}. 

Then C is a (n, N, k, d)-code, where N = N 2 + MNr and d = min{di, d 2 , 2dR}. We write 
C = Cl *Cj{ C -2 for the resulting linkage code and call C the code obtained by linking Ci and C 2 
through Cr. 

The notation C = Ci C 2 has to be used with care because the code C depends on the 
representing sets Mi and M 2 and not only on the codes Ci and C 2 . Thus the notation Mi*Cii -^2 
is actually the accurate one, but we prefer the former because the properties of the linkage that 
we are interested in are associated to the subspace codes Ci, C 2 . The sets Mi, M 2 are mere 
technicalities in our context, and the notation Ci *Cr C 2 will not lead to any confusion. 

However, in order to illustrate the dependence on the SC-representing sets for Ci and C 2 
we will present, after the proof, an example showing that different choices lead in general to 
different distance distributions of the linkage code. 

Proof. The cardinality of C is clear because the three sets Ci are pairwise disjoint. Furthermore, 
it is obvious that ds(Ci) = ds(Cj) for i = 1,2. Moreover, it is clear that each subspace in C 2 
intersects trivially with each subspace in Ci and C 3 . Thus ds(Wi,VV 2 ) = 2k for all Wi G C 2 and 
VV ’2 € Cl U C 3 . 

Next, let U = im([/ | 0) € Ci and V = im(I7' | M) € C 3 . Thus U, U' € Mi and M G Cr\{ 0}. 
Then rk(M) > dR by linearity of the code Cr and so dimker M < k — dR, where kerM = {x G 
F^ I xM = 0}. Let now v gU CiV. Then v = x{U \ 0) = y{U' \ M) for some x, y gV^ and thus 
y G ker M and xU = yU' G im(C/) nim(C/'). Since the maps from UG\V mapping x{U \ 0) to xU 
and y{U' \ M) to y are both injective, this yields 

d\m{U n V) < min{dim(im(C/) PI im(C/')), dim(ker(M))}. 

Now one concludes ds(W, V) > min{(ii, 2dR}, as desired. 

Lastly, let U = im(f7 | M), V = im(f7' | M') G C 3 , and let lA fiV, thus U U' or M M'. 
Let V gU nV. Then v = x{U \ M) = y{U' \ M') for some x, y G F^. Hence 

xU = yU' G im(17) n im(C/') and xM = yM' G im(M) n im(M'). 
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If [/ / U', then dim(W n V) < dim(im(C/) fl im{U')) and ds{U, V) > ds(im(C/), im{U')) > di. 
If [/ = U', then X = y because U has full row rank. Moreover, M — M' G Ci^\{0} and 
thus dim(ker(M — M')) < k — dpt- Now x G ker(M — M') along with the injectivity of the map 
x{U I M) I —X fromZ//nV to ker(M —M') shows that dim(Z//nV) < dim(ker(M —M')) < k — dn, 
and thus ds(Z//, V) > 2dji. This concludes the proof. □ 


The following example shows that different choices of the SC-representing sets for Ci and C 2 
lead to different distance distributions of the linkage code. Since we will not further study the 
distance distribution of linkage codes, we continue to use the notation Ci *Cr C -2 for the linkage. 

Example 2.4. Let {ni,n 2 ,k,q) = (4,4, 2,2) and 


A4i = A42 


{( 


10 1 
0 10 


A 0 0 

Vo 1 0 



and M'i = 



1 

1 



fO 1 0 0\ I 

Vi 1 0 oyj 


and 



0 0 \ 
0 oj 


A 0 0 0\ /O 1 0 0\ 
Vo 1 0 oy’Vo 0 1 oy 




Note that C(A4i) = C(A4{). We find that in the linkage code A4i *Cr-M -2 (see the notation of the 
paragraph after Theorem 12.3p there exist 5 pairs of distinct subspaces with subspace distance 2 
and all other pairs have subspace distance 4, whereas in the linkage code M'l *Cr M 2 only 3 
pairs have subspace distance 2 and all others have subspace distance 4. 


The next two examples illustrate that we can easily construct very large codes of long length 
by suitable linkage. 

Example 2.5. We aim at constructing a constant-dimension code over F 2 of length 13, dimen¬ 
sion k = 3, and distance 4. Let ni = 7 and n 2 = 6 . The largest known codes of dimension 3 
and length 7 (resp. 6 ) with distance 4 have cardinality 329 (resp. 77), see [3l Tables I and II] by 
Braun/Reichelt as well as [28] by Honold et ah, where it is shown that 77 is actually the largest 
possible size for length 6 . We choose these codes for A and C 2 , respectively, and an MRD code Cr 
in with rank distance dpi = 2 and thus cardinality A^j = = 2^^ due to (|2.2I) . The 

resulting linkage code has therefore cardinality 77-|- 2^^ • 329 = 1, 347,661. This is lower than the 
cardinality of the best known code with the same parameters, which is 1, 597, 245, see O App.]. 
But the latter has been found by extended computer search, whereas the linkage code is readily 
available once Ci and A have been found. The linkage code beats the codes that have been found 
by Etzion/Silberstein with the aid of the multilevel construction [TO] where the best such code 
has cardinality 1,192, 587, see [23]. It also beats the modified multilevel construction m by Et¬ 
zion/Silberstein, where the best such code has size 1, 221, 296. This particular construction is a 
refinement of the pending dot construction appearing first in [30] by Trautmann/Rosenthal. The 
following table presents the cardinality of further linkage constructions for g = 2, fe = 3, ds = 4 
and various lengths n. We make use of the best codes found in [3] Table II] for lengths 6 ,..., 9 
and an MRD code Cji C F 3 xn 2 rank distance 2. Thus Npi = 2 ^"^ thanks to (12.2p . We also 
show the largest size obtained via the modified multilevel (MML) construction [111 Thm. 17] 
(which always beats the multilevel construction in [TO]) as well as the largest size known so far. 
For n < 14 the latter has been found by computer search [S] Tables I and II], while for n = 15 
no such search has been conducted yet and linkage with ni = 9, re 2 = 6 results in the largest 
known code. Note that for all n shown in the table, every partition into n = ni + n 2 leads to 
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a linkage code that is larger than the MML construction. This is probably due to the fact that 
the MML construction leads to subspace codes that contain a lifted MRD code. This restriction 
also restricts the size of these codes. 


n 

ni 

n2 

Ni 

N 2 

Linkage 

MML 

Largest Known 

12 

6 

6 

77 

77 

315,469 

305,324 

385,515 

13 

6 

7 

77 

329 

1,261,897 

1,221,296 

1,597, 245 

13 

7 

6 

329 

77 

1,347, 661 

1,221,296 

1,597, 245 

14 

7 

7 

329 

329 

5,390,665 

4, 885,184 

5,996,178 

14 

6 

8 

77 

1,312 

5,047, 584 

4, 885,184 

5,996,178 

14 

8 

6 

1,312 

77 

5,374,029 

4, 885,184 

5,996,178 

15 

8 

7 

1,312 

329 

21,496,137 

19,540, 736 

23,322,701 

15 

7 

8 

329 

1,312 

21,562,656 

19,540, 736 

23,322,701 

15 

6 

9 

77 

5,694 

20,190, 782 

19,540, 736 

23,322,701 

15 

9 

6 

5,694 

77 

23,322,701 

19,540, 736 

23,322,701 


It is worth pointing out that for length n = 13 the largest known cardinality 1,597,245 is 
actually the optimum by the anticode bound, and the existence of a code with that size has 
been established by Braun et al. in [2] via a 2-analogue of a Steiner triple system. 

Example 2.6. Let us consider Theorem 12.31 for the case where Ci is a lifted MRD code and Cr 
is an MRD code. For C 2 we may choose a lifted MRD code or an arbitrary subspace code. Let 
us consider the case where C 2 is a lifted MRD code. Thus, let n = ni -|- n 2 , where rii > 2k 
for i = 1,2, and Ci be the lifting of an MRD code in with rank distance d. Then 

\Ci\ = and ds(Ci) = 2d. Moreover, let Cr C MRD code of rank 

distance d. Thus, = q^^ik-d+i)^ gy Theorem 12.31 the linkage code Ci *Cr C .2 has subspace 
distance 2d and cardinality gi("^ 2 -fc)(fc-d-i-i) _j_ ^{n-k)(k-d+\) ^ Note that the second term is the 
cardinality of a lifted MRD code in with subspace distance 2d. Thus, linkage always results 
in a better code than lifting. In fact, with our choice the code C\ U C 3 in Theorem 12.31 is a 
lifted MRD code and thus the cardinality of the linkage code is clearly larger than that of a 
lifted MRD code. Furthermore we observe that only the first term depends on the partition 
n = ni -|- n 2 , and that the cardinality of C\ *Cr C 2 is largest when n 2 is largest. The following 
table shows the size of the linkage construction for g = 2, fc = 3, ds = 4 and various lengths. In 
each case, Ci is a lifted MRD code of distance 4 and Cr is an MRD code of rank distance 2. Each 
given length n is split into n = ni -|- n 2 such that n 2 is maximal subject to n* > 2 k for i = 1 , 2 . 
In the column denoted by “Linkjargest” we present the cardinality of the linkage code where we 
use the largest known subspace code for C 2 . In the column “LinkMRo” we use a lifted MRD 
code for C 2 . For comparison we also show the size, of a lifted MRD code of length n. 

It should be noted that the linkage codes are smaller than the codes obtained from the MML 
construction; see the previous table. This is explained by the fact that the MML construction 
is a careful design to create additional subspaces without compromising the distance. It may be 
regarded as a replacement of the code C 2 in Theorem 12.31 bv a larger set, where the zero block 
matrix is replaced by suitable matrices. In the column “Extended Lifted MRD” we illustrate 
that our codes are slightly smaller than those constructed in m by SkacheI0, which are also 

^One can show that in our situation the optimal design choice for the parameter hn = in m is hn = 0, 

and therefore the cardinality of the resulting code is as in Section IV.C of [27]. 
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subspace codes containing a lifted MRD code (and are smaller than the MML codes). 


n 

111 

n2 

hinkjargest 

LinkMRD 

Lifted MRD 

Extended Lifted MRD 

12 

6 

6 

262,221 

262,208 

262,144 

266,304 

13 

6 

7 

1,048,905 

1,048,832 

1,048,576 

1,065,216 

14 

6 

8 

4,195,616 

4,195,328 

4,194,304 

4, 260,864 


We will return to these codes in Theorem 14.61 when we investigate decoding. 

The linkage construction can be viewed as a generalization of two specific constructions that 
can be found in the literature. We will discuss the details in the next section, where we turn to 
subspace codes with largest possible distance. 


3 Partial Spreads 


With the aid of Theorem 12.31 we can construct optimal partial spreads for certain cases. Recall 
that a partial spread in F” is a collection of subspaces that pairwise intersect trivially. If all 
subspaces have the same dimension, say k, then this is simply a constant-dimension code of 
dimension k and distance 2k, and we call the code a partial k-spread. It is well known that if k 
divides n, then an optimal partial A:-spread (i.e., a partial A:-spread of maximum cardinality) is a 
k-spread, i.e., the spaces intersect trivially and cover the entire F". In this case a simple counting 
argument shows that the cardinality is (g” — — 1), where F = Fg. Several constructions 

of /c-spreads are known. For later reference we provide the following two options. 

Remark 3.1. Let k divide n. 

(a) |24[ Thm. 11] The orbit of the subfield F^fc in the field F^n under the natural action of the 
group F*n is a fe-spread in F^n. 

(b) |18] Thm. 6, Rem. 8] Let m = n/k and M G GLfc(F) be the companion matrix of a primitive 
polynomial of degree k. Then the set {im(j4i,... ,Am) \ Ai € Fg[M], not all Ai are zero} is 
a A:-spread in F*^. It is called a Desarguesian spread. 


If k does not divide n, then the maximum size of a partial /c-spread in F” is in general not 
known - with one exception which will be considered below in further detail. The following 
result can be found in [ll Thms. 4.1, 4.2] and [71 Thm. 7]; see also | 8 l Thm. 3]. 

Theorem 3.2. Let n (mod k) = c, where c G {0, ...,A: — 1}. Denote the largest possible 
cardinality of a partial k-spread in F” by fi{n, k). Then 

gn _ c 

h^{n,k) > — -- -q^' + l 

q'^ — I 

with equality if c £ {0,1}. Furthermore, if c> 1 then 


! , X ^ q - q 

F[n,k) < —r -- 

g'' — 1 


[ 6 *J — 1 , where 9 


^l + 4gfc(gfc_gc) _ - 2q^ + 1) 

2 
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Constructions of partial fc-spreads and cardinality 


m{n, k) 



+ 1, where n (mod k) = c, 


(3.1) 


were presented in [U Thms. 4.2] as well as |12l Thm. 11] and [191 Thm. 13]. Hence for c G {0,1} 
these partial spreads have maximum possible cardinality. The latter two constructions are 
special cases of our linkage and will be described in our terminology in the following examples. 


Example 3.3. We describe the construction of partial /c-spreads by Etzion/Vardy [121 Thm. 11]. 
Let n >2k and write n = Ik + c, where c € {0,..., A; — 1}. Set ni = k{l — 1) and n 2 = k + c. 
Consider the Fg-vector space F^^i x Fgn 2 . In choose the fc-spread Ci given by the orbit 
of the subfield F^fc under the action of the cyclic group F*ni; see Remark l3.H al. Furthermore, 
let /3 be a primitive element of F^nj and set C 2 = {spanpll, /3,..., Note that this is 

trivially a partial fe-spread in Fg»i 2 of maximal possible cardinality because k > n 2 / 2 . Consider 
the coordinate map w.r.t. the basis {1, / 3 , ..., 1 } of ¥qn 2 , that is, 

n 2 —1 

i=0 

Using the identification ip, the code C 2 simply translates into C 2 = {im(/fc j Okxc)}- Finally, 
in F *^^”2 choose the rank-metric code 


Cr = {{Ik I 0kxc)M^ \j = 0,...,q^^-2}U {0}, 

where M is the companion matrix of the minimal polynomial of /3 over ¥q. Note that the matrix 
{Ik I 0kxc)M^ consist exactly of the rows p{l3^), ■ ■ ■ By Remark 1 2. 11 Cr is an MRD 

code of rank distance k. Identifying Fg^i with F^U the linkage code Ci C 2 is exactly the 
partial /c-spread constructed in [TH Thm. 11]. It has cardinality 1 -|- q^^{q^^ — l)/{q^ — 1), and 
this is m{n, k). 

Example 3.4. Essentially the same construction as in Example 13.31 but with different specifica¬ 
tions of the constituent codes is used by Gorla/Ravagnani in [191 Thm. 13]. Again, \qI n = Ik + c 
and set ni = k{l — 1) and n2 = k + c. Then the code constructed in [191 Thm. 13] is the linkage 
Cl with the following specifications: Ci is a Desarguesian A:-spread in F”i (see Remark LS.ip 

while C2 is the subspace code {im(0fcxc I Ik)} and Cr is an MRD code in F^^”^ 3,3 in Remark 12. II 
with matrix W = {0 \ Ik), thus the nonzero matrices in Cr are the last k rows of the matrices MK 
In addition to the construction, the authors also present a decoding algorithm for their partial 
spreads by making explicit use of the structure of the Desarguesian spread; see [IHl Sec. 5]. In 
contrast, no decoding algorithm is given in m for the partial spreads constructed therein. 

Instead of partitioning n into ni -|-u ,2 with the specific choice of 77-2 = A:-|-c as in the previous 
examples, we may use any other splitting n = Ik + n 2 . This will be summarized in the next 
result where we also address maximality of the partial spread. A partial A:-spread in F*^ is called 
maximal if it is maximal with respect to inclusion, that is, it is not properly contained in any 
other partial A:-spread. The following result shows, among other things, that linking a A:-spread 
and a maximal partial A;-spread through an MRD code leads to a maximal partial A:-spread. 
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Theorem 3.5. Let n = Ik + n 2 , where I > 1 and n 2 > k. Let Ci he a k-spread in and 
C 2 he a partial k-spread in Furthermore, let Cr he a linear MRD code in with rank 

distance k and thus cardinality . Finally, let C = Ci C 2 be the resulting linkage code as in 
Theorem \2.3[ . 

(a) If IC 2 I = m{n 2 , k), then \C\ = m{n, k). 

(b) If C 2 is a maximal partial k-spread then so is C. 

Proof, (a) Theorem 12.31 tells us that \C\ = q^‘^m{lk,k) + m{n 2 ,k). But this is easily seen to 
be m(n, k). 

(b) Let {Wi I W 2 ) G ^kx{m+n 2 ) q£ k and set W := im(lTi | IT 2 ). We have to show 
that there exists a subspace V € C such that W fl V 7 ^ {0}, for then C is a maximal partial 
fc-spread. Assume first that Wi / 0. Then there exists (x, y) G W such that x ^ 0. Since Ci is 
a spread of F**', the vector x is in exactly one subspace of Ci, say im(?7i). Let x = aUi, where 
a G F^\{0}. Since Cr is a linear rank-metric code with rank distance k, we have aM 7 ^ aM' for 
all distinct M, M' G Cr. This shows that the set {aM \ M G Cr} has cardinality \Cr\ = and 
therefore equals F’^^. As a consequence, y = aM for some M G Cr. Hence {x,y) = a{Ui \ M) 
and yy n V 7 ^ {0} for the subspace V = im(C/i | M) G C. Assume now W = im(0 | VL 2 ). Then 
rk(H 2 ) = k and the maximality of C 2 implies WnV 7 ^ {0} for some subspace V = im(0 | U 2 ) G C. 
All of this shows that C is a maximal partial fe-spread in F"". □ 

The following is an immediate consequence of Examples 13.31 and 13.41 because in both cases 
the chosen code C 2 is trivially a maximal partial /c-spread. 

Corollary 3.6. The partial spreads constructed in Thm. 11] and in U9l Thm. 13] are 
maximal. 

Maximality of the partial spreads in [19] has also been established by Gorla/Ravagnani 
in [IHl Prop. 20]. 

In Theorem 13.21 we have seen that the maximum cardinality of a partial fc-spread in F” is 
known whenever n (mod k) G {0,1}. There is one more case where the cardinality is known, 
and that is if g = 2 and k = 3. The following result covers all remainders of n modulo 3. 

Theorem 3.7 ([ 8 l Thm. 5]). Let k = 3 and n > 6 . Let n (mod 3) = c. Then the maximum 
cardinality of a partial 3-spread in Wtf is 

211 _ 2c 


We call a partial 3-spread with this cardinality a maximum partial 3-spread. 

Note that for c G {0,1} the result is simply a special case of Theorem 13.51 whereas for c = 2 
the cardinality m{n, k) in (13.11) is one below the maximum. As a consequence, the constructions 
in [121 Thm. 11] and [191 Thm. 13] are just one subspace short of being maximum. 

The proof of Theorem 13.71 is based on a concrete example for n = 8 and an extension 
construction for n > 8. It makes use of a result in [U Lem. 4], which establishes a partition 
of Fg into subspaces of two distinct dimensions. Below we will provide an alternative extension, 
where we will also make use of the maximum 3-spread in F®. 
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Example 3.8 ([ 8 l Ex. 2]). There exists a partial 3-spread in F® with cardinality 34. Hence the 
spread is maximum. It has been found by computer search and is explicitly given in [ 8 ]. 

Now we can provide a simple construction of maximum partial 3-spreads in F 2 for any n > 
10. Note that, due to the previous example and earlier discussions, a maximum partial 3-spread 
in F 2 is available for the values n G { 6 , 7, 8 ,9}. 

Corollary 3.9. Let n > 10 and write n = 31 + n 2 for some I > 1 and n 2 G {6, 7, 8}. Choose a 
3-spread Ci in F^^ and a maximum partial 3-spread C 2 in F^^. Finally, let Cr he an MRD code 
with rank distanee 3 in Then Ci *Cr C 2 Is a maximum partial 3-spread in F^. 


Proof. The resulting code is certainly a partial spread. Let n (mod 3) = c, thus n 2 (mod 3) = c. 
By Theorems 13.71 and 12.31 the cardinality of Ci *Cr C -2 is given by 


23 ' - 1 2*^2 _ 2 = 

2 +^-c = 

7 7 



and this is the maximum value due to Theorem Ml 


□ 


4 Decoding of Linkage Codes 

In this section we turn to decodability of the linkage codes from Theorem 12.31 Of course, one 
aims at reducing decoding of Ci *CrC2 to decoding of the smaller codes Ci, C 2 , Cr. We will show 
first that this strategy does not work if one utilizes the rank metric for the code Cr. Instead 
one has to employ the subspace distance for all codes involved. We will show that if we use 
suitable MRD codes and liftings thereof, then decoding can indeed be reduced to decoding of 
the constituent codes. Since lifted Gabidulin codes can be efficiently decoded, as proven by Silva 
et al. [26], this leads to an efficient decoding algorithm for a particular instance of linkage codes. 
The following terminology is standard. 

Definition 4.1. Let C be a subspace code in F"" with subspace distance d. A subspace V C F*^ 
is called decodahle w.r.t. C if there exists a subspace U ^C such that ds(7/, V) < 

Since the subspace distance is a metric on the set of all subspaces in F*^ (see [221 Lem. 1]), 
a decodable subspace has a unique closest codeword in C. 

The following simple fact will be useful later. 

Remark 4.2. Let C be a constant-dimension code in F” with dimension k and subspace dis¬ 
tance d. Suppose V C F” is a decodable iL-dimensional subspace and 7/ G C is the unique closest 
codeword. Then 

dim(7/ n V) > K/2 and dim(W + V) < A: + 77/2. 

The first inequality follows from ds(7f, V) < (d—1)/2 < k, see (|2.ip . which then reads as K + k — 
2dim(7fnV) < k. The second inequality is obtained by using dim(7f-|-V) = A:-|-i7 —dim(7/n V). 

We start with an example illustrating that the rank distance of the code Cr cannot be used 
in the natural way for decoding the linkage code Ci C 2 . 

Throughout this section we call a matrix V G F^^*^ a matrix representation of the sub¬ 
space V C F"' if V = im(M) := {xM \ x G F^}. We explicitly allow dim(V) < K. 
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Example 4.3. Let ni = 77-2 = 8, /c = 4 and F = F 2 . Define 


Cl = C 2 = {im(/ I I), im(I | 0)}. 


where I and 0 are the identity and the zero matrix in F^^^, 


respectively. Moreover, let 


Cr = { 04 x 8 , (C I 0), (M I 0), (/ + M I 0)}, where M 


/O 1 0 0\ 
0 0 10 
0 0 0 1 
\1 1 0 0 / 


Then ds(Ci) = ds(C 2 ) = 8 and cIr := dR(CR) = 4. Thus the linkage code C = Ci *Cjf C 2 has 
length 16, subspace distance d = 8 and cardinality N = 10. Consider the received word 


V = im 


1100110001100000 \ 

oiiooiioooiiooooy' 


Then ds(V, U) = 2 for the codeword U = im([/i | U 2 ) G C, where 


Ui = {I\ I) and C /2 = (M I 0). 

In particular, V is decodable. Note that V ^lA (i.e., only erasures occurred during the trans¬ 
mission). One can check straightforwardly that there exists no matrix representation {Vi \ V 2 ) 
of V in F^^^® such that dR(V 2 , C/ 2 ) < = |. Even worse, for all matrix representations 

(Li I V 2 ) G F^^^® of V for which the matrix V 2 bas a unique closest matrix in Cr with re¬ 
spect to the rank distance, this unique closest matrix is the zero matrix and therefore does 
not lead to the correct decoding U. The fact that the “obvious decoding” does not work may 
be explained by the fact that the subspaces represented by the nonzero matrices in Cr, i.e., 
im(I I 0), im(M | 0), im(/ -|- M | 0), all coincide. This causes the rank-metric code Cr to be of 
little help with decoding. 


The example can be generalized. We introduce the following notation. Define the projections 

TTi : F’^i+^2 ^ F^^i, (a, b) ^ a and 7r2 : F’^i+’^2 ^ ^ (41) 

For a subspace M C F"" we define Ai = 7ri{A). Thus, if A = im(74i | ^42), then Ai = im(^j) for 

i = 1,2. 

Proposition 4.4. Let C be as in Theorem 1^.31 and assume d > Ar + 2. Then there exists a 
subspace U = im(C/i | C/2) G C and a received word V C F"" such that 

(a) ds{U,V) < (that is, V is decodable), 

(b) V(LU (hence only erasures occurred during transmission), 

(c) for any V 2 G F^^”^ im(V2) = 7r2(V) we have 

rk(E2 - U 2 ) > 

In other words, it is not possible to decode V by making use of the rank metric for the code Cr. 
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Proof. Let d > dn + 2. Since d < di < 2k, we have r := k — and thus r > 1 because 

it is an integer. We construct now subspaces U and V as stated in the proposition. First choose 
a subspace U = ini([/i | U 2 ) G C with {Ui \ U 2 ) € such that rk(C/ 2 ) > k — r = [%■]. By 

definition of C such an element does indeed exist. Next, there exists a matrix X € GLfc(F) such 
that 

and where U 22 G j^as rank k — r and and im(f72i) Ll im(C/ 22 ) = {0}. Put 


V = im{Uu \U2i). 

Then dim(V) = r because the rows of {Ui \ U 2 ) are linearly independent. Moreover, V ^lA and 
therefore 

ds(V,^) = A: + r - 2 r = 

This establishes (a) and (b). For (c) consider now all matrices V 2 in F^^”^ whose row space is 
7 r 2 (V). These matrices can be written as 

V 2 := ^ F^^*” is any matrix of rank r. 

The matrix X~^ does not change the row space, and we include it only to simplify the next 
step. Indeed, for each such matrix V 2 we have 

rk(P 2 - U 2 ) = rk(XP 2 - XU 2 ) = rk - ^ 21 ^ > rk(M 2 t/ 2 i - t/ 22 ). 

The rightmost matrix has full row rank, k — r. Indeed, suppose u{M 2 U 2 i — U 22 ) = 0. Then 
UM2U21 = UU22 € im(t/2i) n im([/22)- Since this intersection is trivial, we obtain UU22 = 0, 
which in turn implies u = 0. All of this shows that rk ([/2 ~ P 2 ) > for all matrix 

representations of 7 r 2 (V). □ 

The last observation suggests to modify the linkage construction by simply replacing the 
rank-metric code Cr by matrix representations of a subspace code. This results in a code that is 
decodable if its constituent codes are decodable. But since these codes are considerably smaller 
than the original linkage codes, we will not follow that path. 

Instead, we will show now how to decode linkage codes Ci *Cr C 2 for the case where Ci and Cr 
are (lifted) MRD codes. We need the following lemma. 

Lemma 4.5. Let C be the code from Theorem \2.d[ and V C F*^ he a deeodable K -dimensional 
subspace. LetU be the closest eodeword, thus ds(W, V) < Then 

lAi = 0 dimVi < K/2. 

Proof. Assume U.\ 7 ^ 0. Then rk([/i) = A; by dehnition of C. Thus dim(W) = k = dim(Wi), 

and the map tti\i( is injective, where tti is the projection from (|4.1I1 . We compute 

dim(Z/A n V) = dim( 7 ri(W H V)) < dim( 7 ri(Z/A) n 7 ri(V)) < dim( 7 ri(V)) = dimVi < K/2, 
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which is a contradiction to Remark 14.21 Thus = 0. 

Let Ui = 0, thus U = im(0 | U 2 ) and rk(? 72 ) = khy definition of C. Write V = im(i/i | V 2 ) 
for some (Vi | V 2 ) € With the aid of Remark 14.21 we obtain 

rkCLi) + A; < rk ( ° = dim{U + V) <k + —. 

V 2 J 2 

Hence dimVi = rk(Vi) < K/2. □ 

Note that the implication of the last lemma is in general not true for U 2 and V 2 

because the matrices in Cr may not have rank k. 

Now we are in the position to discuss decoding of the linkage codes from Theorem 12.31 We 
consider the following situation which is a special case of the general linkage construction. It 
may also be regarded as an extension of the codes considered in Example 12.61 

Theorem 4.6. For i = 1, 2 let rii > k and let A4j C be a linear MRD code with 

rank distance d, thus \Mi\ = Moreover, let C F^^"i and M. 4 ^ C 

SC-representing sets of constant-dimension codes with subspace distance 2d. Consider the code 
C = C' U C" U C", where 


C = {im(4 I Ml I M 2 ) I Ml E Ml, M 2 E M 2 }, 

C" = {im(0fcxfc I M I Ofcxna) I M E M^}, 

C'" = {im(0A;xfc I Ofcxni I M) I M E A 44 }. 

Then C is an (n, N, k, 2d)-code, where n = fc + ni + n 2 and N = q{'n-i+n2){k-d+i) _|_ |_y\ 42 | _|_ \M 4 \. 

Proof. This is a simple application of Theorem l2.3l : C = Ci *^^^2 with the codes Ci = {im(/fc)}, 
C 2 = {im(M I 0) I M E M 3 } U {im(0 | M) | M E M 4 }, and Cr = {(Mi | Ms) | M* E Mi}. □ 

Note that if A44 represents a lifted MRD code, then C' U C'" = Ci *m 2 ^ 2 , where Ci = 
{im(I I Ml) I Ml E Ml}, C 2 = {im(M) | M E M 4 }. Hence the code is of the form as discussed 
in Example 12.61 and in Theorem 14.61 we improve upon the codes in that example by the size 
of M 3 . In the column “LinkMRo” of the table in Example 12.61 we listed, for a specific choice 
of parameters, the largest codes of the form C' U C" above. As we saw already, the largest size 
is attained when ns is largest subject to n = A: + ni + ns with ni + A: > 2A: and ns > 2 k (now 
ni + A; takes the role of ni from that example), thus for ni = k and ns = n — 2 k. But in that 
case I Ada I = 1 and thus we improve upon the codes in that table by exactly one subspace. 

We now turn to decoding of the codes in Theorem 14.61 As we show next this can be reduced 
to decoding of the constituent codes. 

Theorem 4.7. Consider the setting of Theorem \4-(^ To ease notation we set no := k. For 

i = 1,2 define the lifted MRD codes Ci := {im(/fc | M) | M E Ad*}. Furthermore, set C 3 := 

^(Mfa) and C 4 := C(M 4 ). Then, if C\,..., C 4 are decodable then so is the linkage code C. More 
precisely, let V = im(Eo | Li | V 2 ) C F”, Vi E be a K-dimensional received word such that 

ds(V,C) < {2d — l)/2. Then exactly one of the following situations occurs. 

(a) rk(Eo | Li) < 7^/2. In this case the unigue closest codeword in C is in C" and given by U = 
im(0 I 0 I M), where M E Ad 4 is the unigue matrix such that ds(im(M), im(E 2 )) < (2d—1)/2. 


13 


(b) rk(Vb I V2) < i^/2. In this case the unique closest codeword in C is in C and given by U = 
im(0 I M I 0), where M G is the unique matrix such that ds(im(M), im(Vi)) < (2d—1)/2. 
(^cJrk(Vo) > K/2. In this case the unique closest codeword in C is in C and given by U = 
im(I I Ml I M2), where Mi G M.i are the unique matrices such that ds(im(/ | Mj), im(Vo | Vi)) < 
(2d - l)/2 fori = 1,2. 


Proof. First of all, the uniqueness of the matrices M and Mi in (a) - (c) is guaranteed since the 
subspace codes C{M. 3 ), and the lifted MRD codes Ci all have subspace distance 2d. 

Let us denote hy U = im(t/o I Ui \ U 2 ) the unique codeword in C closest to V. We will use 
the notation Ui = im([/j) and Vj = im(Vi) for i = 0 , 1 , 2 . 

First we show that at most one of the 3 cases above can occur. Clearly, if rk(Vo) > iL/2, 
then neither (a) nor (b) can occur. Let now rk(Vo) < iL/2 and assume r := rk(Vo | Vi) < iL/2. 
We have to show that rk(Vo | V 2 ) > K/2. After suitable row operations we may assume 

(Vol U I U, = (''»■ K-). 

where the first block row has r rows. Then rk(Vo | Vi | V 2 ) = K implies that rk(V 22 ) > Ar/2. This 
implies rk(Vo | V 2 ) > K/2, as desired. Using symmetry, all of this shows that if rk(Vo) < K/2, 
then at most one of the cases (a) or (b) can occur. 

Next we show that exactly one of the cases (a) - (c) occurs. To do so, it suffices to show 
that if rk(Vb) < K/2 then (a) or (b) must occur. We know that C = Ci C 2 with Ci, C 2 , Cr 
as in the proof of Theorem 14.61 Therefore, Lemma 14.51 along with rk(Vo) < K/2 implies Uq = 0. 
Thus, U £ C" U C". Let us assume U G C, say U = im(0 | M | 0). Along with rk(M) = k we 
derive 

rk(Uo I U 2 ) + < rk ^ = dim(L/ + V) < fc + K/2, (4.2) 

where the last inequality is due to Remark 14.21 As a consequence, rk(Vo | U 2 ) < K/2. Simi¬ 
larly, U G C" implies rk(Vo | Ui) < Ar/2. 

Now we turn to decoding for each of the three cases. 

(a) Let rk(Vb | Ui) < K/2. Suppose the closest codeword U is in C, say U = im(0 | M | 0). 
Then (14.2p shows that rk(Vo | V 2 ) < K/2. But this means that also case (b) occurs, a contradic¬ 
tion. Hence U £ C Li C". But this code is a linkage code. Indeed, 

C'uC"' = Ci *c«C4, 

where Ci = {im(/fc|Mi) | M\ £ A 4 i}, C4 = C(A44), and Cr = M.2. Hence Lemma 14.51 
implies (I/q | Ui) = (0 | 0 ), thus U = im (0 | 0 | M) G C" with M £ AI4. In particular, rk(M) = k. 
Consider the projection 7r2 of onto F"'^. Then {'it2)\u is injective and thus dim(WnV) < 

dim(ZY2 n V2). This implies ds(V2, U2) < K + k — 2dim(V2 n U2) < ds(V, U) < { 2 d — l)/ 2 . 
Thus, decoding V2 to its closest codeword in C{Mi) results in U2. Using its unique matrix 
representation M £ M.4, i.e., U2 = im(M), we arrive at the correct decoding U = im (0 | 0 | M) 
of the received space V. 

(b) The case rk(Vb | U 2 ) < K/2 is analogous. 

(c) Let rk(Vo) > K/2. Then Lemma [4.51 applied to Ci *Cr C 2 (see proof of Theorem 14.61) 

implies (C/q) 0. Thus U £ C. In particular, we may assume C/q = Ik- For i = 1,2 let V) := 
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im(Vb I Vi) and U- := im(/fc | Ui). Then W' € C* for z = 1,2, where Ci = {im(/fc | M) | M € 
is the lifting of the MRD code Mi for z = 1, 2. Consider the projections 

^ (00,01,02) ^ (ocO^) 

Then (z/;*) \u is injective and thus d\m.{U n V) < dim(Z^' n V'). As in (a) this implies ds(Z^(, V') < 
ds(Z^, V) < i2d — l)/2 for z = 1,2. Hence V' can be uniquely decoded w.r.t. Ci and the closest 
codeword is given by U[. Using the unique matrix representations (/|[/j), Ui G Mi, of the 
spaces Ui, we arrive at the correct decoding of V. □ 

We summarize the result in the following algorithm. 


Algorithm 4.8: Decoding algorithm for the codes in Theorem 14.61 
Data: a decodable A-dimensional subspace V = im(Vb | Ui | V2) with (Vq | Ui | U2) € 
Result: the unique U G C' U C" U C" such that ds(V,W) < 
if rk(Vb I Ui) < ^ then 

decode im(V2) in C{Mi) to im(C/2); 
return U = im (0 | 0 | C/ 2 ), 
else 

if rk(Vb I U 2 ) < -y then 

decode im(Vi) in ClMs) to im(C/i); 
return U = im (0 | C/i | 0 ). 
else 

decode im(Vo | Ui) in Ci to im(/fc | C/i); 
decode im(Vo | V2) in C2 to im(/fc | U2); 
return U = im(/fc | Ui | U2). 


One should observe that in the last case of the algorithm, the two decoding steps can be 
performed in parallel. A similar, but not identical, form of parallelizing decoding is also used 
for the spread codes in [19]; recall Example 13.41 for the relation to our linkage codes. 

Clearly, the construction in Theorem 14.61 and its decoding can easily be generalized to more 
than 3 blocks. 

Remark 4.9. A very efficient decoding is obtained when we use Gabidulin codes for A4i and M 2 
and lifted Gabidulin codes for M 3 and A44 (thus rzj > k). In this case, all codes relevant for 
decoding in the previous proof are lifted Gabidulin codes, and the decoding algorithm derived 
by Silva et al. [26] may be employed. If zzj » k, then even better efficiency is obtained by using 
direct products of Gabidulin codes as the MRD codes and the lifting of such a code for M 3 
and Ma', see [261 Sec. VI.E]. In fact, our code C in Theorem 14.61 (or rather its generalization 
to more than 3 blocks) is of the form proposed in |26j and with the above we have shown how 
to enlarge the code without compromising its properties. In this sense, our results put the 
considerations in [261 Sec. VI.E] in a broader context. 
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